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We summarize the results of two recent searches for flavor-changing neutral current, lepton- 
flavor violating, and lepton-number violating decays of D + , Df , and D° mesons (and their 
antiparticles) into modes containing muons and electrons. Using data from Fermilab charm 
hadroproduction experiment E791, we examined D + and nil and Kit decay modes and 
the D° dilepton decay modes containing either £ + £~ , a p°, K , or <f> vector meson, or a non- 
resonant 7T7T, Kit, or KK pair of pseudoscalar mesons. No evidence for any of these decays was 
found. Therefore, we presented branching-fraction upper limits at 90% confidence level for 
the 51 decay modes examined. Twenty-six of these modes had no previously reported limits, 
and eighteen of the remainder were reported with significant improvements over previously 
published results. 



1 Introduction 

The E791 Collaboration has previously reported limits on rare and forbidden dilepton decays of 
charged charm mesonsE™. Such measurements probe the SU(2) xU(l) Standard Model of electroweak 
interactions in search of new mediators and couplings @a. Here we summarize the results of two 
related analyses. First i we examined the nil and Kit decay modes of D + and Df and the £ + £~ 
decay modes of D°. Then we extend the methodology to 27 dilepton decay modes of the D° meson 

1 containing either resonant V£ + £~~ decays, where V is a p°, K*°, or <fi, and non-resonant hxhzll 
decays, where hi is either a tt or a K. The leptons were either muons or electrons. Charge-conjugate 
modes are implied. The modes are lepton flavor- violating (e.g., D + — ► 7r + /^ + e~), or lepton number- 
violating (e.g., Df — * 7r _ /i + /i + ), or flavor-changing neutral current decays tag., D° — » K e + e~). 
Box diagrams can mimic FCNC decays, but only at the 1(T 10 to 1(T 9 level Hfl. Long range effects 
through resonant modes (e.g., D° — > K*° p°, p° — > e + e~) can occur at the 10 -6 leveli'0. Numerous 
experiments have studied rare decays of charge -1/3 strange quarks. Charge 2/3 charm quarks are 
interesting because they may exhibit a different couplings. 

The data come from measurements made with the Fermilab E791 spectrometer H. A total of 

2 x 10 10 events were taken with a loose transverse energy requirement. These events were produced by 
a 500 GeV/c n~ beam interacting in a fixed target consisting of five thin, well-separated foils. Track 
and vertex information came from "hits" in 23 silicon microstrip planes and 45 wire chamber planes. 
This information and the bending provided by two dipole magnets were used for momentum analysis 
of charged particles. Kaon identification was carried out by two multi-cell Cerenkov counters 1121 that 
provided n/K separation in the momentum range 6 — 60 GeV/c. We required that the momentum- 
dependent light yield in the Cerenkov counters be consistent for kaon-candidate tracks, except for 
those in decays with <j> — > K + K~ , where the narrow mass window for the <f> decay provided sufficient 
kaon identification (ID). 

Electron ID was based on transverse shower shape plus matching wire chamber tracks to shower 
positions and energies in an electromagnetic calorimeter 0. The electron ID efficiency varied from 
62% below 9 GeV/c to 45% above 20 GeV/c. The probability to misidentify a pion as an electron was 
~ 0.8%, independent of pion momentum. 

Muon ID was obtained from two planes of scintillation counters. The first plane (5.5 m x 3.0 m) 
of 15 counters measured the horizontal position while the second plane (3.0 m x 2.2 m) of 16 counters 
measured the vertical position. There were about 15 interaction lengths of shielding upstream of the 
counters to filter out hadrons. Data from D + — ► K decays 113 were used to choose selection 

criteria for muon candidates. Timing information from the smaller set of muon scintillation counters 
was used to improve the horizontal position resolution. Counter efficiencies, measured using muons 
originating from the primary target, were found to be (99± 1)% for the smaller counters and (69±3)% 
for the larger counters. The probability of misidentifying a pion as a muon decreased with increasing 
momentum, from about 6% at 8 GeV/c to 1.3% above 20 GeV/c. 

Events with evidence of well-separated production (primary) and decay (secondary) vertices were 
selected to separate charm candidates from background. Secondary vertices were required to be 
separated from the primary vertex by greater than 20<t l for D + decays and greater than 12 a L for D° 
and Df decays, where a L is the calculated resolution of the measured longitudinal separation. Also, 
the secondary vertex had to be separated from the closest material in the target foils by greater than 
5 a' L , where a' is the uncertainty in this separation. The vector sum of the momenta from secondary 
vertex tracks was required to pass within 40 pm. of the primary vertex in the plane perpendicular 
to the beam. The net momentum of the charm candidate transverse to the line connecting the 
production and decay vertices had to be less than 300 MeV/c for D candidates, less than 250 MeV/c 
for Df candidates, and less than 200 MeV/c for D + candidates. Finally, decay track candidates were 
required to pass approximately 10 times closer to the secondary vertex than to the primary vertex. 
These selection criteria and kaon identification requirements were the same for both the search mode 
and for its normalization signal (discussed below). 



To determine our selection criteria, we used a "blind" analysis technique. Before the selection 
criteria were finalized, all events having masses within a window AM5 around the mass of the D° 
were "masked" so that the presence or absence of any potential signal candidates would not bias our 
choice of selection criteria. All criteria were then chosen by studying events generated by a Monte 
Carlo (MC) simulation programEj and background events, outside the signal windows, from real data. 
The criteria were chosen to maximize the ratio Nmc/V^b, where Nmc and Ng are the numbers of 
MC and background events, respectively, after all selection criteria were applied. The data within the 
signal windows were unmasked only after this optimization. We used asymmetric windows for the 
decay modes containing electrons to allow for the bremsstrahlung low-energy tail. The signal windows 
were: 1.83 < M(D°) < 1.90 GeV/c 2 for /x/x and 1.76 < M(D°) < 1.90 GeV/c 2 for ee and fie modes. 

The upper limit for each branching fraction Bx was calculated using the following formula: 



D = A^fNorm x B where fNorm = ( , } 

AWm ex e x fx 



Nx is the 90% confidence level (CL) upper limit on the number of decays for the rare or forbidden decay 
mode X and B^f OTm is the normalization mode branching fraction obtained from the Particle Data 
Group 111. £Norm and ex are the detection efficiencies while f$Sm an( ^ /i° are the fractions of Monte 
Carlo events that were reconstructed and passed the final selection criteria, for the normalization and 
decay modes, respectively. 

The 90% CL upper limits Nx are calculated using the method of Feldman and Cousins 113 to 
account for background, and then corrected for systematic errors by the method of Cousins and 
Highland @. In these methods, the numbers of signal events are determined by simple counting, not 
by a fit. Upper limits are determined using the number of candidate events observed and expected 
number of background events within the signal region. (See Refs. iH for a more detailed discussion of 
backgrounds.) 



2 The D + -» Ml, Df -> hit and D° -> i+£~ Analysis 

We normalized the sensitivity of our search to topologically similar Cabibbo-favored decays. For the 
D+ decays we used 24010 ± 166 D+ K-ir+ir + ; for Df decays we used 782 ± 30 Df -> 4>ir + ; and 
for D° decays we used 25210 ± 179 D° -> K~ir + events. The widths of our normalization modes were 
10.5 MeV/c 2 for D + , 9.5 MeV/c 2 for Df, and 12 MeV/c 2 for D°. The results are shown in Table | 
and compared with previous results in Figure 111. 



Table 1: E791 90% confidence level branching fraction upper limits for D + -> Ml, Df -> Ml and D° -> l + T. 



Mode E791 Limit 


Mode E791 Limit 


Mode E791 Limit 


D+ TT + n + ii- 1.5 x 10~ 5 
D + TT-fl + fI + 1.7 x 10~ 5 
D+ _» K+fx+fj,- 4.4 x 10~ 5 


D + -> 7r + e + e" 5.2 x 10~ 5 
D + -> vr-e+e + 9.6 x 10~ 5 
L>+ -> A"+e+e- 2.0 x 10" 4 


D + -► vr+^e^ 3.4 x 10~ 5 
D+ -► vr^/i + e + 5.0 x 10~ 5 
£)+ _> K+^e^ 6.8 x 10~ 5 


Df -► K+n+fi- 1.4 x 10~ 4 
Df -► K-/J.+H+ 1.8 x 10~ 4 
Df -» vr+^ + fT 1.4 x 10" 4 
Df vr-^+^+ 8.2 x 10" 5 


L>+ -> K + e + e- 1.6 x 10~ 3 
L>+ -> A"- e +e + 6.3 x 10~ 4 
L>+ -> vr+e+e" 2.7 x 10~ 4 
L>+ -> 7r-e + e+ 6.9 x 10" 4 


D+ K + fi ± e T 6.3 x 10" 4 
L>+ K~fi + e + 6.8 x lO" 4 
D+ -» vr+^e^ 6.1 x 10~ 4 
Df vr-/i+e + 7.3 x 10" 4 


D° -► n + fi~ 5.2 x 10" 6 


D° e+e" 6.2 x 10~ 6 


L>° -> /r^ 8.1 x 10~ 6 



3 The D° VW" and £>° -> Analysis 

There were a few minor differences between this analysis and our previous analysis as discussed above. 



First, we examined resonant modes, where the mass ranges used were: 



< 150 MeV/c 2 



m K - 



< 55 MeV/c 2 , and \rrix+K- ~ m <t>\ < 10 MeV/c 2 . We normalized the sensitivity of 
each search to similar hadronic 3-body (resonant) or 4-body (non-resonant) decays. One exception is 
the case of D° -» p°£ ± £+ where we normalize to nonresonant D — * tt+tt tt+tt because no published 
branching fraction exists for D° — > p°tt+tt- . Table || lists the normalization mode used for each signal 
mode and the fitted numbers of normalization data events (iVNorm)- 



Table 2: Normalization modes used for D 



o 



Vl+T and D° -» hhll. 



Decay Mode Norm. Mode A^orm 


Decay Mode Norm. Mode Nj<f OTtn 


D° -> p ^^ L>° -> Tr+vr-vr+TT- 2049±53 

L>° -> L>° -> 07T+7T- 113±19 

I? -> ETvr^ L>° -> K-tt + tt-tt + 11550±113 


D° -> iTV^ £)0 _^ ^*o^+^- 5451±72 
L>° -> ttttU D° -> it + it-it + it- 2049±53 
D° -> L>° -> K+i^-vr+vr- 406±41 



The final results are shown in Table || and compared with previous results in Figure 



Table 3: E791 90% confidence 



level branching fraction upper limits for D° -> Vl+T and D° 



hhll. 



Mode E791 Limit 


Mode E791 Limit 


Mode E791 Limit 


D° -> TT+TT- Ll+Ll' 3.0 X 10" 5 

L>° -> K-TT+Li+fi- 3.6 x 10~ 4 
L>° -> K+K- Ll+Ll- 3.3 x 10~ 5 


D° -» vr+vr-e+e- 3.7 x 10" 4 
L>° K-vr+e+e" 3.9 x 10" 4 
D° K+i^-e+e" 3.2 x 10" 4 


D° -» TT+TT' Li+e+ 1.5 X 10~ 5 

D° K-ir+f^e* 5.5 x 10" 4 
L>° K+K-fi+e* 1.8 x 10" 4 


L>° -> p°ll+li- 2.2 x 10~ 5 
D° K*° Li+Li- 2.4 x 10" 5 
D° -> <MV~ 3.1 x 10~ 5 


L>° -» p°e+e- 1.2 x 10" 4 
D ^K*°e+e- 4.7 x HT 5 
L>° 4>e+e- 5.9 x 10~ 5 


L>° -» p°/i ± e =F 6.6 x 10~ 5 
D° -► K*°Li ± e+ 8.3 x 10~ 5 
L>° 4.7 x 10~ 5 


L>° -> TT-TT-LL+LL+ 2.9 X 10~ 5 
L>° -> K-TT-Ll+fL+ 3.9 x 10~ 4 

L>° -> K-R-LL+LL+ 9.4 x 10" 5 


L>° -» vr-7r-e + e+ 1.1 x 10" 4 
L>° K-7r-e+e + 2.1 x 10" 4 
L>° K-i^-e + e+ 1.5 x 10" 4 


D° -» TT~ TT~ LL+ €+ 7.9 X 10~ 5 

D° K-7r-^+e + 2.2 x 10" 4 
L>° K-R-LL+e+ 5.7 x 10~ 5 




Figure 1: Comparison of the 90% CL upper-limit 
branching fractions from E791 data (dark circles) 
with existing limits (open diamonds) from the 1998 
PDG 14 . 



Figure 2: Comparison of the 90% CL upper-limit 
branching fractions from E791 data (dark circles) 
with existing limits (open diamonds) from the 2000 
PDG 14 . 



4 Conclusion 

We used a "blind" analysis of data from Fermilab experiment E791 to obtain upper limits on the 
dilepton branching fractions for 51 flavor-changing neutral current, lepton-number violating, and 
lepton- family violating decays of D + , Df, and D° mesons. No evidence for any of these 2, 3 and 
4-body decays was found. Therefore, we presented upper limits on the branching fractions at the 
90% confidence level. Eighteen limits represented significant improvements over previously published 
results. Twenty-six of the remaining modes had no previously reported limits. 
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